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Abstract 

We analyze Monte Carlo simulation and series-expansion data for 
the susceptibility of the three-state Potts model in the critical re- 
gion. The amplitudes of the susceptibility on the high- and the low- 
temperature sides of the critical point as extracted from the Monte 
Carlo data are in good agreement with those obtained from the se- 
ries expansions and their (universal) ratio compares quite well with a 
recent quantum field theory prediction by Delfino and Cardy. 



1 Introduction 



The universal thermodynamic behaviour of a system in the vicinity of a 
critical point is characterized by a set of critical exponents and by universal 
combinations of critical amplitudes [0]. While the critical exponents are 
generally well studied, there are still few theoretical results on the universal 
combinations of critical amplitudes. Recently, some progress was achieved by 
Delfino and Cardy [Q] (this reference shall be denoted as article I throughout 
our paper) for the two-dimensional g-state Potts model and some universal 
amplitude ratios were computed for g = 2, 3 and 4 (see also for the most 
recent review on other results.) 

It is commonly accepted that for a typical spin model, for instance the 
Ising model, there are only two independent length scales in the criti- 
cal region and thus there must be four universal relations |^ among the six 
following critical amplitudes: the amplitudes of specific heat in the ordered 
phase A- and in the symmetric phase the analogous amplitudes of the 
magnetic susceptibility r_ and r+; the amplitude of the magnetization when 
approaching the critical temperature from below, -B; and the amplitude of the 
dependence of the magnetization on the magnetic field at the transition tem- 
perature, D. In the case of the two-dimensional Potts model with g > 2, it 
is possible to define also a "transverse susceptibility" in the low temperature 
phase 0. In terms of its critical amplitude F^, a fifth universal ratio Tt/T- 
can be defined which is determined only by the behaviour in the LT phase. 

The values of the ratio r_|_/r„ of the susceptibility critical amplitudes, 
were calculated in paper I with the results 37.699, 13.848 and 4.013 for q = 
2, 3 and 4, respectively. The first value coincides with the well known exact 
result in the four digits presented. In Ref. ^ Delfino, Barkema and Cardy 
performed a Monte Carlo (MC) test of the other predictions of paper I and 
found results for r+/r_ not consistent with their expectations in the q = 3 
case, while in the g = 4 case their results were inconclusive. Actually the 
analysis of the data in the g = 4 case is somewhat difficult due to the expected 
logarithmic corrections & 0, ffO] to the power-like critical behaviour and the 



results [Tl[] of a multi-parameter fit are still controversial (see, for instance 
Ref. [|^). On the other hand their analytical calculations for Tt/T- are in 
very good agreement with the numerical results reported in Ref. 0. 

Here we present an analysis of the susceptibility of the g = 3 Potts model 
by a MC simulation supported by extrapolations of the presently available 
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LT and HT series expansions. The results of both procedure are completely 
consistent in the critical region window. The simplest extrapolation of the 
series expansions by Pade approximants yields the estimate r+/r_ = 14.2(3), 
while a fit of the MC data leads to r+/r_ = 14 ± 1. Thus the agreement 
with the Delfino and Cardy prediction for g = 3 is very good. Since also 
their estimate of the susceptibility ratio for g = 2 is correct, our result gives 
us further confidence that even the prediction for the g = 4 case might be 
correct. However, further numerical study of this case would still be welcome. 

The text is organised as follows. In section ^ we define the model and 
present a few details of the series calculations and of the MC simulations. 
The data fitting procedure is discussed in the section |^ and the results are 
summarised and commented in section |^. 

2 Model and Definitions 

The Hamiltonian of the Potts model is 



where Sj is a site variable taking q integer values between and g — 1, and h 
is an external magnetic field which stabilizes the state 0. 

The susceptibility is defined in terms of the free energy per site h), 
where (3 = 1/T is the inverse temperature variable, starting with a finite 
system of sites and then taking the thermodynamic limit: 



H= ^(l-5,^,,J + /,5](l-5,^,o) 



(1) 



<ij> i 



f{P,h)= hm -l^\og{ZN{P,h)) 



(2) 



with a partition function Zn defined according to 



(3) 



conf 



The susceptibility at inverse temperature (3 is thus 



(3 dh^ 



h=0 



(4) 
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2.1 Low-temperature and high-temperature expansions 



The first few LT expansion coefficients can be very simply obtained by clas- 
sifying tlie spin configurations with respect to their energy and multiphcity. 
In the LT variable z = we have 



X={q- + 2{q - 1)^^ + 2{q - l){q - 2)z' + ... (5) 



It is by far less trivial to derive the very long LT expansions through z' 
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tabulated in Ref. [1^ for zero field in the q = 2,3 and 4 cases. 



On the other hand, the presently available HT expansions for the suscep- 



tibility are still of moderate length. For general g, they reach [T^ the order 
10 in the HT variable 



l + iq-l)z' 

The first few terms of the expansion for x are 

X = l + Av + 12v^ + 36v^ + (76 + 12q)v^ + ... 

The HT series is however sufficient to compute accurately x provided that 
one does not get too close to the critical point Pc = log(l + ^g). It should 
be noted that the HT expansion tabulated in Ref. is conventionally 
normalized to 1 at infinite temperature and therefore it should be multiplied 
by the factor to get the correct amplitude. 

We recall that the critical exponent is given by the exact formula 

Ttt^ - 8/i7r + 

IW = 



with = 2 cos /i. 

In terms of the HT series and of the known values of 7 and f3c, we can 
form, on the HT side of the critical point, the "HT effective amplitude" 
r+(r) = r^x(r). Here r = 1 — /3//3c is the reduced inverse temperature. 
Similarly, we can construct a "LT effective amplitude" r_(r) = (— r)^x(r), 
on the LT side r < of the critical point. The critical amplitudes r+ and 
r_ are computed by extrapolating to the critical point the corresponding 
effective amplitudes. Since the HT series is not very long, the extrapolation 
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can only be performed, in the most naive way, by Pade approximants, which 
cannot allow for the singular corrections to scaling. Therefore we should not 
expect an accuracy of the amplitude r_|_ better than a few percent. The 
LT series is much longer, but we have preferred to use Pade approximants 
also in this case. We have used the highest order available approximants to 
extrapolate the effective amplitudes, namely the [5/5] approximant in the 
HT region, leading to r+ = 0.180(3) and the [22/22] approximant in the 
LT region, which gives r_ = 0.0127(1). In both cases we have checked the 
accuracy of the results by comparing these approximants with the lower order 
ones which are nearby in the Pade table. From this study we can conclude 
that r+/r_ = 14.2(3). Since the largest contribution to the error comes 
from the determination of r_|_, in order to improve the result, significantly 
longer HT series should be computed and the analysis should be performed 
by differential approximants, which can allow for the corrections to scaling. 



2.2 Monte Carlo simulations 

We adopt the Wolff algorithm [jl^ for studying square lattices of linear size 



L with periodic boundary conditions. Starting from a typical ordered state, 
we let the system relax in 10^ steps measured by the number of flipped Wolff 
clusters. The averages are computed over 10^ steps. Our random numbers are 
produced by an exclusive- XOR combination of two shift-register generators 
with the taps (9689,471) and (4423,1393), which are known to be safe 
for the Wolff algorithm. 

The maximal system size L = 200 we used is rather moderate for the 
current computer standards, but is quite sufficient for our purpose, as we 
will see. 

During the simulations we have evaluated an order parameter M{t) de- 
fined as follows 

- 1 

M = ^ (6) 

q-1 ^ ^ 

where Nm is the number of sites i with Sj = m at the time t of the sim- 



ulation ||T6|, and m G [0...{q — 1)] is the value of the majority of the spins 
at that time. This is like using the modulus of the magnetization in a MC 
simulations of Ising model. 



4 



Thus, the susceptibihty in the LT phase is given by the fluctuations of 
the majority of the spins 

ksTx- = ^(< iV^ > - < iV™ >') (7) 
while in the HT phase x is given by fluctuations in all q states, 

kBTx+ = E -^(< Nl>-<N, >2), (8) 

where Ni is the number of sites i with the spin in the state /i. Properly 
allowing for the finite-size effects, this definition of the susceptibility gives, 
in both phases, an extremely good consistency with the series expansion data. 



3 Data analysis 

In this section we shall argue that system sizes such that L x L ^ 100 x 100 
sites are sufficiently large to reproduce the susceptibility behaviour of the 
infinite system in the range of the "critical temperature window" . Then we 
shall present a comparison of MC and series data. Finally, we shall discuss 
the fit to the data and comment on the quality of the results. 



3.1 Critical region window 



The critical region window |r| G [Tfs,Tcorr] is usually characterized 
as an interval between two (inverse reduced) temperatures tjs and Tcorr in 
which the finite system of size L exhibits the critical behaviour of the infinite 
system. For instance, in the case of the susceptibility we should have 



XOcr±|r|-^ (9) 

where r± are the critical amplitudes in the HT phase (+) and in the LT phase 
(-), respectively, and 7 is the critical exponent. The "left-hand" boundary of 
the window Tfg is the rather well defined temperature at which the correlation 
length ^ cx r^*^ is of the order of the lattice size L, so that tjs ^ L~^^'^. 
The "right-hand" limiting value of the critical window Tcorr is related to the 
corrections to scaling and could be defined as the temperature at which the 
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deviation from the critical behaviour due to the corrections to scahng 
reaches the level of e.g. 1%. Therefore the temperature Tcorr could be well 
identified knowing exactly the corrections to scaling as in the case of the 
Ising model [1^, |18[. In the case of insufficient theoretical knowledge of 
the amplitudes of the correction terms, as is the case of the 3-state Potts 
model, we can try to identify the "right-hand" temperature boundary Tcorr 
by plotting the MC data for the previously defined "effective amplitude" . 

The MC results are shown in the Figure |l| together with the series expan- 
sion data which are denoted by thick solid lines. The series data compare 
well with the MC data also for size L < 80 when |r| > 0.01. The corrections 
to scaling are not small in the interval of the reduced temperature accessible 
by the Monte Carlo simulations. 

3.2 Fit to the data 

We have fitted the data taken in the critical window of the LT and HT phases 
to the following expression 

X = r±|r|-^(l + a±|r|^)+D± (10) 



where 7 = 13/9 and A = ?/t,2/i/t,i = 2/3 are known exactly []I9|, The value 
of A is supported by the series expansion analysis of Ref . ||2^ , by a finite-size 
analysis |]2T]] of the 3-state Potts quantum chain, and by a finite-size analysis 



of the transfer matrix results of Ref. [^]. The constant "background" terms 



D± are known to be important even for the Ising model |]T3, especially so 
in the LT phase. We have kept a single correction-to-scaling term a±|r|^ in 
order to avoid introducing too many fitting parameters. However, we have 
also checked the stability under inclusion of further terms, e.g., (1 + a±|r|^ -|- 
^±|t| + c±|rp^) and found that T± varies only within the accuracy of the fit. 

The results of the fit to the susceptibility data are shown in Table |I] for 
lattices sizes L = 20, 40, 60, 80, 100 and 200. The right-hand boundary of the 
critical window was chosen as |r| = |(T — Tc)/T\ = 0.2 and the left-hand 
boundary according to the lattice size as |r| > l/L^^". 

For comparison, we have also fitted in the same way the series expansion 
data, taking as critical window the interval |r| G [0.01 : 0.2] and assuming 
that the series expansion data become increasingly accurate for larger values 
of |r|. The error bars refiect the accuracy of the fit and do not allow for 
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Figure 1: Finite-size behaviour of the effective amphtude of the suscep- 
tibiUty data in the critical region window. The data for lattice sizes 
L = 20,40,60,80,100 and 200 are denoted by squares, left-triangles, di- 
amonds, down-triangles, up-triangles and circles, respectively. The high- 
temperature data correspond to the'^upper part and the low-temperature 
data to the lower part of the figure. The data generated from the high- and 
low-temperature series expansions are shown by thick solid lines. 



Table 1: Results of the fit to the MC susceptibility data in the critical region 
window. For comparison, we have also reported in the last line the results of 
a similar fit to the series expansion (SE) data. 



L 


r_ 


a_ 


D_ 


r+ 


a+ 


D+ 


r+/r_ 


20 


0.01321(3) 


-1.61(1) 


0.0087(5) 


0.2036(6) 


-0.67(2) 


0.37(1) 


15.41(9) 


40 


0.01260(2) 


-1.500(9) 


0.0071(3) 


0.1966(3) 


-0.64(1) 


0.448(8) 


15.60(5) 


60 


0.01245(4) 


-1.375(1) 


0.0026(1) 


0.1642(3) 


0.90(2) 


-0.228(6) 


13.19(7) 


80 


0.01251(1) 


-1.393(7) 


0.0026(3) 


0.1815(2) 


0.14(1) 


0.008(5) 


14.51(3) 


100 


0.01250(1) 


-1.41(6) 


0.0037(3) 


0.1728(16) 


0.49(9) 


-0.11(5) 


13.82(16) 


200 


0.01273(1) 


-1.509(6) 


0.0070(2) 


0.1741(2) 


0.46(1) 


-0.17(1) 


13.68(3) 


SE 


0.012774(3) 


-1.517(8) 


0.0070(2) 


0.1783(7) 


0.24(2) 


0.005(6) 


13.96(7) 



systematic deviations due to the fact that the critical window is far from the 
asymptotic limit when large corrections to scaling are present. The fit to 
the MC data shows stability and consistency with the series expansion data, 
listed in the last line of Table |1|. 

Very conservatively, we can conclude from our MC data and Fade approx- 
imation of series expansion that the ratio of the susceptibility amplitudes for 
the 3-state Fotts model is F+/F_ = 14 ± 1. Thus, our results are com- 
pletely consistent with the value F+/F_ = 13.848 calculated by Delfino and 
Cardy §. 

An additional check of the results is obtained by studying the ratio 
r*^'^^)(r) of the HT and the LT effective amplitudes of the susceptibility, as 
computed from series expansions r('^'^)(|r|) = F+(r)/F_(— r). We can expect 
the following behaviour of this ratio 

r^^'WA) = ^ (l + ailr^ + a,\r\ + a,\T\'^ + ...) (11) 

as r — > 0+. 

A three-parameter (^, ai, and 02) fit of r'^^^\\T\) in the temperature window 
[0.01 : 0.1] gives ^ = 14.2(5). 

In the case of the Ising model (g = 2 Fotts model), a similar fit of the 
same ratio r'^^^\\T\) to the form 
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r(^^)(|r|) = ^ (1 + + a2|r| log |r|) (12) 

leads to a value of the critical amplitude ratio ^ = 39±2 which is consistent 
with the exact value 37.69365.... The large error bars are due to the fact that 
we have used series expansions of the same length as for the 3-state Potts 
model in order to test under similar conditions the accuracy of the method. 
Of course, much better results could obtained fully using the much longer 
series expansions which are available for the 2d Ising model ||2^ . 



4 Summary of the results and conclusions 

The values of the suceptibility critical amplitude ratio for the Potts model 
with q = 2,3 and 4 were calculated by Delfino and Cardy using the two- 
kink approximation of the exact scattering theory for the Potts model [P^]. 
The value 37.699 thus obtained for the ratio in the q = 2 case (Ising model) 
agrees well with the exactly known value. However the same authors and 
Barkema were unable to confirm numerically ^ their theoretical results: 
13.848 for the g = 3 Potts model and 4.013 for the 4-state Potts model. The 
discussion of the latter case is beyond present paper. However by analysing 
our MC data and the existing series expansions, we find that the critical 
amplitude ratio for the 3-state Potts model can be very safely identified with 
the estimate 14 ± 1, quite consistently with the prediction by Delfino and 
Cardy. 

What is the main difference between the analysis of Ref. and that 
presented here? First, we calculate the amplitudes separately in both the 
LT and the HT phases by fitting the temperature behaviour of the suscep- 
tibilities. It is also important that the value of susceptibility was computed 
by not less than 10^ Wolff MonteCarlo steps at each value of temperature. 
Indeed the fit to the susceptibility becomes unstable for smaller statistics. 
Next, we have computed the amplitude ratio as a function of temperature 
defined in the same way in both phases. This is not the case for the analy- 
sis in Ref. [0], where the corresponding temperatures in the two phases are 
shifted by a factor proportional to the ratio of the correlation lengths. 

Since two out of the three r_|_/r_ ratio values computed by Delfino and 
Cardy agree well, either with the known exact result for the Ising case, or 
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with our MC and series expansion data for g = 3, little doubt remains, in 
our opinion, that also their prediction for the 4-state Potts model may be 
correct. We wish to quote here a MC analysis of the 4-state Potts model 
by Caselle et al. where an estimate consistent with the Delfino and Cardy 
prediction is obtained. However Delfino et al. did not found this analysis 
completely satisfactory and therefore the susceptibility ratio prediction for 
g = 4 is still waiting for further numerical verifications. 



5 Acknowledgements 

LNS is grateful to the Theoretical Physics group of the Milano-Bicocca Uni- 
versity and to the Statistical Physics group of the University Henri Poincare 
Nancy-I for their kind hospitality. Financial support from the twin research 
program between the Landau Institute and the Ecole Normale Superieure 
de Paris as well as financial support from the CARIPLO Foundation and 
Landau Network — Centro Volta, and Russian Foundation for Basic Research 
under project 99-07-18412 are also acknowledged. 



10 



References 

[1] V. Privman, P.C Hohenberg, A. Aharony, in Phase Transitions and 
Critical Phenomena, Vol. 14, edited by C. Domb and J.L. Lebowitz 
(Academic, New York, 1991). 



[2] 
[3] 
[4] 
[5] 
[6] 

[7] 

[8] 
[9] 

[10 

[11 
[12 

[13 
[14 
[15 
[16 
[17 
[18 



G. Delfino and J.L. Cardy, Nud. Phys. B 519 [FS] (1998) 551. 



G. Mussardo, |tiep-th/0010T64 . 



D. Stauffer, M. Ferer, and M. Wortis, Phys. Rev. Lett. 29 (1972) 345. 

A. Aharony and P.C. Hohenberg, Phys. Rev. B 13 (1976) 3081. 

T.T. Wu, B.M. McCoy, C.A. Tracy and E. Barouch, Phys. Rev. B 13 
(1976) 316. 

G. Delfino, G.T. Barkema, and J.L. Cardy, Nucl. Phys. B 565 (2000) 
521. 

N. Nauenberg and D.J. Scalapino, Phys. Rev. Lett. 44 (1980) 837. 

J. L. Cardy, N. Nauenberg and D.J. Scalapino, Phys. Rev B 22 (1980) 
2560. 

J. Salas and A. Sokal, J. Stat. Phys. 88 (1997) 567. 

M. Caselle, R. Tateo, S. Vinti, Nucl. Phys. B 562 (1999) 549. 

K.M. Briggs, I.G. Enting, and A.J. Guttmann, J. Phys. A 27 (1994) 
1503. 

G. Shreider and J.D. Reger, J. Phys. A. 27 (1994) 1071. 
U. Wolff, Phys. Rev. Lett. 62 (1989) 361. 
L.N. Shchur, Comp. Phys. Comm., 121-122 (1999) 83. 
K. Binder, J. Stat. Phys, 24 (1981) 69. 

A.L. Talapov and L.N. Shchur, J. Phys.: Cond. Matt., 6 (1994) 8295. 
L.N. Shchur and O.A. Vasilyev, to be published in Phys. Rev. E 

11 



[19] B. Nienhuis, J. Phys. A, 15 (1982) 199. 

[20] J. Adler and V. Privman, J. Phys. A, 15 (1982) L417. 

[21] G. von Gehlen, V. Rittenberg, and T. Vescan, J. Phys. A, 20 (1987) 
2577. 

[22] S.L.A. de Queiroz, J. Phys. A, 33 (2000) 721. 

[23] W. P. Orrick, B. Nickel, A. J. Guttmann, and J. H. H. Perk, Phys. Rev. 
Lett. 86 (2001) 4120; J. Stat. Phys. 102 (2001) 795. 

[24] L. Chim and A.B. Zamolodchikov, Int. J. Mod. Phys. A, 7 (1992) 5317. 



12 



